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Abstract Breakup reactions are one of the main tools for the study of exotic nu- 
clei, and in particular of their continuum. In order to get valuable information from 
measurements, a precise reaction model coupled to a fair description of the pro- 
jectile is needed. We assume that the projectile initially possesses a cluster struc- 
ture, which is revealed by the dissociation process. This structure is described by 
a few-body Hamiltonian involving effective forces between the clusters. Within 
this assumption, we review various reaction models. In semiclassical models, the 
projectile-target relative motion is described by a classical trajectory and the reac- 
tion properties are deduced by solving a time-dependent Schrodinger equation. We 
then describe the principle and variants of the eikonal approximation: the dynami- 
cal eikonal approximation, the standard eikonal approximation, and a corrected ver- 
sion avoiding Coulomb divergence. Finally, we present the continuum-discretized 
coupled-channel method (CDCC), in which the Schrodinger equation is solved with 
the projectile continuum approximated by square-integrable states. These models 
are first illustrated by applications to two-cluster projectiles for studies of nuclei far 
from stability and of reactions useful in astrophysics. Recent extensions to three- 
cluster projectiles, like two-neutron halo nuclei, are then presented and discussed. 
We end this review with some views of the future in breakup-reaction theory. 
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1 Introduction 

The advent of radioactive ion beams has opened a new era in nuclear physics by pro- 
viding the possibiHty to study nuclei far from stability. In particular the availability 
of these beams favoured the discovery of halo nuclei [1 j. Due to the very short life- 
time of exotic nuclei, this study cannot be performed through usual spectroscopic 
techniques and one must resort to indirect methods. Breakup is one of these meth- 
ods. In this reaction, the projectile under analysis dissociates into more elementary 
components through its interaction with a target. Many such experiments have been 
performed with the hope to probe exotic nuclear structures far from stability |2i|3j|. 

In order to get valuable information from breakup measurements, one must have 
not only a fair description of the projectile, but also an accurate reaction model. 
At present, a fully microscopic description of the reaction is computationally un- 
feasible. Simplifying assumptions are necessary. First, we will discuss only elastic 
breakup, i.e. a dissociation process leaving the target unchanged in its ground state. 
Other channels are simulated through the use of optical potentials. Second, we as- 
sume a cluster structure for the projectile. The projectile ground state is assumed to 
be a bound state of the clusters appearing during the breakup reaction. The bound 
and continuum states of the projectile are thus described by a few-body Hamilto- 
nian involving effective forces between the constituent clusters. Theoretical reaction 
models are therefore based on this cluster description of the projectile and effective 
cluster-cluster and cluster-target interactions. 

Even within these simplifying model assumptions, a direct resolution of the re- 
sulting many-body Schrodinger equation is still not possible in most cases. In this 
article, we thus review various approximations that have been developed up to now. 

We begin with the models based on the semiclassical approximation H in which 
the projectile-target relative motion is described by a classical trajectory. This ap- 
proximation is valid at high energies. It leads to the resolution of a time-dependent 
Schrodinger equation. In a primary version, the time-dependent equation was solved 
at the first order of the perturbation theory l?). Then, as computers became more 
powerful, it could be solved numerically |l5]|6]|7l|8]|9][l0l. We present both versions 
indicating their respective advantages and drawbacks. 

We then describe the eikonal approximation flT\ and its variants. The principle 
is to calculate the deviations from a plane-wave motion which are assumed to be 
weak at high energy. By comparison with the semiclassical model, it is possible to 
derive the dynamical eikonal approximation (DEA) that combines the advantages 
of both models lfT2l ITSll . The standard eikonal approximation is obtained by mak- 
ing the additional adiabatic or sudden approximation, which neglects the excitation 
energies of the projectile. With this stronger simplifying assumption, the final state 
only differs from the initial bound state by a phase factor This approach is mostly 
used to model reactions on light targets at intermediate and high energies. Its draw- 
back is that the Coulomb interaction leads to a divergence of breakup cross sections 
at forward angles. This problem can be solved using a first-order correction of the 
Coulomb treatment within the eikonal treatment. A satisfactory approximation of 
the DEA can then be derived lfT4l[T5]| : the Coulomb-corrected eikonal approxima- 
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tion (CCE), which remains valid for breakup on heavy targets. It reproduces most of 
the results of the DEA, although its computational time is significantly lower |fT6l 
which is important for the study of the breakup of three-cluster projectiles. 

Finally, we present the continuum-discretized coupled-channel method (CDCC) 
ifTTlfTSl . in which the full projectile-target Schrodinger equation is solved approxi- 
mately, by representing the continuum of the projectile with square-integrable states. 
This model leads to the numerical resolution of coupled-channel equations, and is 
suited for low- as high-energy reactions. 

All the aforementioned models have been developed initially for two-body pro- 
jectiles. However, the physics of three-cluster systems, like two-neutron halo nuclei, 
is the focus of many experimental studies and must also be investigated with these 
models. We review here the various efforts that have been made in the past few years 
to extend breakup models to three-cluster projectiles 1 1 9 20 2 1] . 

In Sec. 12] we specify the general theoretical framework within which the projec- 
tile is described. The semiclassical model and approximate resolutions of the time- 
dependent Schrodinger equation are described in Sec. [3] Sec. |4]presents the eikonal 
approximation as well as the related DEA and CCE models. Next, in Sec. |5] the 
CDCC method is developed. In Sec.|6] we review applications of breakup reactions 
to two-body projectiles. In particular, we emphasize the use of breakup to study nu- 
clei far from stability and as an indirect way to infer cross sections of reactions of 
astrophysical interest. Sec.|7]details the recent efforts made to extend various reac- 
tion models to three-body projectiles. We end this review by presenting some views 
of the future in breakup-reaction theory. 



2 Projectile and reaction models 

We consider the reaction of a projectile P of mass mp and charge Zpe impinging 
on a target T of mass mr and charge Zpe. The projectile is assumed to exhibit 
a structure made of clusters with masses nij and charges Z,e (mp = and 
Zp = Y-i^i)- Its internal properties are described by a Hamiltonian Hq, depending 
on a set of — 1 internal coordinates collectively represented by notation ^ . With 
the aim of preserving the generality of the presentation of the reaction models, we 
do not specify here the expression of Hq. Details are given in Sees. |6] and |2l where 
applications for the breakup of two- and three-body projectiles are presented. 

The states of the projectile are thus described by the eigenstates of Hq. For total 
angular momentum J and projection M, they are defined by 

i/o<?>^^(£,^)=£0r(£,^), (2.1) 

where E is the energy in the projectile centre-of-mass (cm.) rest frame with respect 
to the dissociation threshold into clusters. Index T symbolically represents the set 
of all additional quantum numbers that depend on the projectile structure, like spins 
and relative orbital momenta of the clusters. Its precise definition depends on the 
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number of clusters and on the model selected when defining Hq. We assume these 
numbers to be discrete, though some may be continuous in some representations 
when there are more than two clusters. To simplify the notation, the parity k of the 
eigenstates of Hq is understood. In the following, any sum over / implicitly includes 
a sum over parity. 

The negative-energy solutions of Eq. (12.1b correspond to the bound states of the 
projectile. They are normed to unity. The positive-energy states describe the broken- 
up projectile with full account of the interactions between the clusters. They are 
orthogonal and normed according to {(j)-^!^ {E' ,^)\^:('^ {E,^)) = d{E — ii')^^^'. To 
describe final states when evaluating breakup cross sections, we also consider the 
incoming scattering states 0-' They correspond to positive-energy states of Hq 

describing the clusters moving away from each other in the projectile cm. frame 
with specific asymptotic momenta and spin projections. These momenta are not 
independent, since the sum of the asymptotic kinetic energies of the clusters is the 
positive energy E. However, within that condition, their directions and, if > 2, 
their norms can vary. By k^, we symbolically denote these directions and wave 
numbers, as well as the projections of the spins of the clusters. These incoming 
scattering states are thus solutions of the Schrodinger equation 



//oC'(£,^)=£C (^,^). (2.2) 



They can be expanded into a linear combination of the eigenstates ^^J^ of Eq. i2.\\ 
with the same energy as 

ir\E,^)^Y^a-'f{k^)^i^{E,^), (2.3) 

^ JMt 

where the coefficients al^ depend on the projectile structure. These scattering states 
are normed following (^i,"' ^ {E,E,)) = 5{E -E')8{k^- k'^). 

The interactions between the projectile constituents and the target are usually 
simulated by optical potentials chosen in the literature or obtained by a folding 
procedure. Within this framework the description of the reaction reduces to the res- 
olution of an (A^+ l)-body Schrodinger equation 
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L.+Ho + Vpt{^,R) 



W{^,R)=EjW{^,R), (2.4) 



where R = {R,Qr) = {R, 9r,(Pr) is the coordinate of the projectile centre of mass 
relative to the target, P is the corresponding momentum, ji = mpmr/imp + mT) 
is the projectile-target reduced mass, and Ej is the total energy in the projectile- 
target cm. frame. The projectile-target interaction VpT is expressed as the sum of 
the optical potentials (including Coulomb) that simulate the interactions between 
the projectile constituents and the target. 
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VpT{^,R)^Y.^iTiRiT), (2.5) 

1=1 

where R/j- is the relative coordinate of the projectile cluster / with respect to the 
target. 

The projectile being initially bound in the state ^tq*'" of negative energy £0, we 
look for solutions of Eq. ( 12.4b with an incoming part behaving asymptotically as 

f (^,R)^^ e'{^^+'"°W«-2)l}0^»^«(£o,<^), (2.6) 

where Z is the component of R in the incident-beam direction. The wavenumber K 
of the projectile-target relative motion is related to the total energy Ej by 

£t = -^+£o. (2.7) 
2m 

The P-T Sommerfeld parameter is defined as 

ri^ZpZTe-lTiv, (2.8) 

where v = hK/ ji is the initial P-T relative velocity. 

A first idea that may come to mind is to solve Eq. i2.4\ exactly, e.g., within the 
Faddeev framework or its extensions. However, the infinite range of the Coulomb 
interaction between the projectile and the target renders the standard equations 
ill-defined. Only recently significant progress has been made. For example, in 
Refs. II22I |23l , this problem is tackled by using an appropriate screening of the 
Coulomb force. This technique has been used to successfully describe the elastic 
scattering and breakup of the deuteron on various targets. However, it has long been 
limited to light targets (see Ref. 1241 for a recent extension to a heavier target). To 
obtain a model that is valid for all types of target, one must still resort to approxima- 
tions in the resolution of Eq. (I2.4l i. These approximations are made in the treatment 
of the projectile-target relative motion, like in the semiclassical (Sec.[3]l or eikonal 
(Sec. lU approximations, or by using a discretized continuum, like in the CDCC 
method (Sec.|5]i. 



3 Semiclassical approximation 

3.1 Time-dependent Schrodinger equation 

The semiclassical approximation relies on the hypothesis that the projectile-target 
relative motion can be efficiently described by a classical trajectory R(f) ID. It is 
thus valid when the de Broglie wavelength is small with respect to the impact pa- 
rameter b characterizing the trajectory, Kb ^ 1, i.e. when the energy is large enough. 
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Along that trajectory, the projectile experiences a time-dependent potential V that 
simulates the Coulomb and nuclear fields of the target. The internal structure of the 
projectile, on the contrary, is described quantum-mechanically by the Hamiltonian 
Hq. This semiclassical approximation leads to the resolution of the time-dependent 
equation 

m^W{^Xt) = [//o + V(^,r)]f (^,b,f). (3.1) 

The time-dependent potential is obtained from the difference between the projectile- 
target interaction Vpj ( 12.51 ) and the potential Vtiaj that defines the classical trajectory 

y(^,r) yfrK,R(f)] - Vtraj[-R(0]- (3.2) 

The potential Viraj acts as a P-T scattering potential that bends the trajectory, but 
does not affect the projectile internal structure. Its interest lies in the fact that V 
decreases faster than VpT- Its effect amounts to changing the phase of the wave 
function. Usually it is chosen to be the Coulomb potential between the projectile 
centre of mass and the target, but it may include a nuclear component. At sufficiently 
high energy, the trajectory is often approximated by a straight line. 

For each impact parameter b, Eq. (13.1b has to be solved with the initial condition 
that the projectile is in its ground state, 

f(^'')(^,b,r)^^^</.^»''«(£o,^). (3.3) 

For each trajectory, the time-dependent wave function 'f'(*'o) must be calculated for 
the different possible values of Mq. 



3.2 Cross sections 

From the output of the resolution of Eq. (13. Il l, the probability of being in a definite 
state of the projectile can be obtained by projecting the final wave function onto the 
corresponding eigenstate of Hq. One can for example compute the elastic scattering 
probability 

This probability depends only on the norm of the impact parameter b because the 
time-dependent wave function f ('^o) depends on the orientation of b, i.e. on the 
azimuthal angle <Pr, only through a phase that cancels out in the calculation of Pel- 
From this probability, the cross section for the elastic scattering in dkection £i is 
obtained as 
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where b{Q) is given by the classical relation between the scattering angle and the 
impact parameter derived from potential Viay The factor do^^^ /dQ, is the elastic 
scattering cross section obtained from Vtraj- In most cases da^^^ /dO. is generated 
from the Coulomb interaction and is thus the P-T Rutherford cross section. 

Likewise, a general breakup probability density can be computed by projecting 
the final wave function onto the ingoing scattering states of Hq, 

■Ph^^) - -^i:i(</.(;'(£,^)|f ('^"'(^-b.r ^ +oo))\\ (3.6) 
dk^dE /Jo + i 

After integration and summation over , the breakup probabihty per unit energy 
reads 

= yrnl^ ^ i(<?..^^(£,^)|f (^°H^,b,f ^ +^)}\\ (3.7) 

dE 2yo + lMoyMT 

Similarly to Eq. ( 13.5b . a differential cross section for the breakup of the projectile is 
given by 

da,.. _^<"^^u[,(^)]. (3.8) 



dEdn dn dE 

The breakup cross section can then be obtained by summing the breakup probability 
over all impact parameters 

da,. 2;rr%(^)M/,. (3.9) 



dE Jo dE 

Because of the trajectory hypothesis of the semiclassical approximation, the im- 
pact parameter b is a classical variable. Therefore, no interference between the dif- 
ferent trajectories can appear This is the major disadvantage of that technique since 
quantal interferences can play a significant role in reactions, in particular in those 
which are nuclear dominated. 



3.3 Resolution at the first order of the perturbation theory 

In the early years of the semiclassical approximations, Eq. ( 13.1b was solved at the 
first order of the perturbation theory [4]. This technique, due to Alder and Winther, 
was applied to analyze the first Coulomb-breakup experiments of halo nuclei [25 1|. 

The time-dependent wave function f ('^o' is expanded upon the basis of eigen- 
states of Hq in Eq. (12.1b . At the first order of the perturbation theory, the resulting 
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equation is solved by considering that V is small. With the initial condition ( 13.31 1. 
the wave function at first order is given by | 



1 



1 + ^ / eT,"^^' V{^,t')e-T,"^^' dt' 
ill J-oo 



Following Eq. ( 13.6b . the general breakup probability density reads 



./.^'^"(Eo,^). (3.10) 



dktdE 



i-2 



2Jo- 



1^ 

^ Mo 



'iE,^)\V{^,t)\^if%Eo,^))dt 



(3.11) 



where (O — (E — Eo)/h. The breakup probabihty per unit energy reads 



dE 



-,-2 



{b) = 



2Jo + 



^ MqJMt 



e"''{rr''iE,^)\V{^,t)\^if%Eo,^))dt 



(3.12) 



With Eq. ( 13. lot , exact expressions can be calculated when considering a purely 
Coulomb P-T interaction for straight-line trajectories in the far-field approximation 
ll27ll . i.e. by assuming that the charge densities of the projectile and target do not 
overlap during the collision. One obtains 

2^^— ^L^T-T^^M(«'^)('/'r(£,^)Kf (^)l0'f°(£o,^)>,(3.13) 

where are the electric multipoles operators of rank X, and I^^ are time inte- 

grals (see, e.g., Eq. (13) of Ref. Il28l ) that can be evaluated analytically as ll27l 



2A + 1 1 



(Ob\ 



(3.14) 



where K„ is a modified Bessel function ||29l . 

If only the dominant dipole term El of the interaction is considered, the breakup 
probabihty (13.12b reads Ii30il 



dP^^ 

dE ^ ' 



167r ( Zje 
Tiv 



(ob 

V 



0)b 

v 



dB{El) 
dE 



(3.15) 



The last factor is the dipole strength function per energy unit 1301 , 
dB{El) 



dE 



2Jo 
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= L L \{(^fiE,^)\^^\^Mif''iEo,m'- (3.16) 

Since modified Bessel functions decrease exponentially, the asymptotic behaviour 
of dP^^/dE for b — > oo is proportional to e\p{—2(ob/v). 

In the case of a purely Coulomb P-T interaction, the first order of the perturbation 
theory exhibits many appealing aspects. First, it can be solved analytically. Second, 
the dynamics part (Ii^) and structure part (matrix elements of ^^^) are separated in 
the expression of the breakup amplitudes ( 13.131 ). This first-order approximation has 
therefore often been used to analyze Coulomb-breakup experiments by assuming 
pure El breakup (see, e.g., Ref. 1251 ). However, as will be seen later, higher-order 
and nuclear-interaction effects are usually not negligible, and a proper analysis of 
experimental data requires a more sophisticated approximation. 



3.4 Numerical resolution 

The time-dependent Schrodinger equation can also be solved numerically. Various 
groups have developed algorithms for that purpose H] |6] |7] |8] |9l [lOl [31] [32l . They 
make use of an approximation of the evolution operator U applied iteratively to the 
initial bound state wave function following the scheme 

f('^o)(^,b,f + 4f) = f/(f + 4f,f)'F(^o)(^,b,0. (3.17) 

Although higher-order algorithms exist (see, e.g., Ref. ||33]| ). all practical calcula- 
tions are performed with second-order approximations of U . Various expressions 
of this approximation exist, depending mainly on the way of representing the time- 
dependent projectile wave function. However they are in general similar to iTOl 

U{t + At) = e-'w.n^-t+A,)^-ifH,^-,%v[^,)^^^^p^_ (3 

With this expression, the time-dependent potential can be treated separately from 
the time-independent Hamiltonian Hq, which greatly simplifies the calculation of 
the time evolution when the wave functions are discretized on a mesh ifTOl . 

The significant advantage of this technique over the first order of perturbation 
is that it naturally includes higher-order effects. Moreover, the nuclear interaction 
between the projectile and the target can be easily added in the numerical scheme 
ll34l . However, the dynamical and structure evolutions being now more deeply en- 
tangled, the analysis of the numerical resolution of the Schrodinger equation is less 
straightforward than its first-order approximation. The numerical technique is also 
much more time-consuming than the perturbation one. The first order of the pertur- 
bation theory therefore remains a useful tool to qualitatively analyze calculations of 
Coulomb-dominated reactions performed with more elaborate models. Moreover, 
as will be seen in Sec. 14.41 it can be used to correct the erroneous treatment of the 
Coulomb interaction within the eikonal description of breakup reactions. 
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Fig. [U illustrates the numerical resolution of the time-dependent Schrodinger 
equation for the Coulomb breakup of ' ^ Be on lead at 68 Me V/nucleon ifTOl . It shows 
the breakup cross section as a function of the relative energy E between the '"^Be 
core and the halo neutron after dissociation. The full line corresponds to the calcu- 
lation with both Coulomb and nuclear P-T interactions. The dashed line is the result 
for a purely Coulomb potential, in which the nuclear interaction is simulated by an 
impact parameter cutoff at bmm = 13 fm. A calculation performed with an impact 
parameter cutoff at b^im = 30 fm simulating a forward-angle cut is plotted as a dot- 
ted line. The experimental data from Ref. 1251 are multiplied by a factor of 0.85 as 
suggested in Ref. Il35l after a remeasurement. 




0.5 1 1.5 2 2.5 3 

E (McV) 



Fig. 1 Breakup cross section of " Be on Pb at 68 Me V/nucleon as a function of the relative energy 
E between the '"Be core and the neutron. Calculations are performed within the semiclassical 
approximation with or without nuclear interaction 1101 . Experimental data 125) are scaled by 0.85 
1351 . Reprinted figure with permission from Ref. 1101 . Copyright (2003) by the American Physical 
Society. 

This example shows the validity of the semiclassical approximation to describe 
breakup observables in the projectile cm. frame for collisions at intermediate en- 
ergies. It also confirms that for heavy targets the reaction is strongly dominated by 
the Coulomb interaction. The inclusion of optical potentials to simulate the nuclear 
P-T interactions indeed only slightly increases the breakup cross section at large 
energy. This shows that Coulomb-breakup calculations are not very sensitive to the 
uncertainty related to the choice of the optical potentials. Nevertheless, since opti- 
cal potentials can be very easily included in the numerical resolution of the time- 
dependent Schrodinger equation, they should be used so as to avoid the imprecise 
impact-parameter cutoff necessary in purely Coulomb calculations. 
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4 Eikonal approximations 

4.1 Dynamical eikonal approximation 

Let us now turn to a purely quantal treatment providing approximate solutions of the 
Schrodinger equation V.A\ . At sufficiently high energy, the projectile is only slightly 
deflected by the target. The dominant dependence of the (A^ + 1 ) -body wave function 
f on the projectile-target coordinate R is therefore in the plane wave contributing 
to the incident relative motion ( |2.6l l. The main idea of the eikonal approximation 
is to factorize that plane wave out of the wave function to define a new function f 
whose variation with R is expected to be small lim[30l[36ll 

f(^,R)=e*^'?(^,R). (4.1) 

With factorization ( 14. Il l and energy conservation ( I2.7l l. the Schrodinger equation 
( 12.41 ) becomes 



— + vPz + //o - £o + (<^ , R) 

2/i 



f(^,R)=0, (4.2) 



where the relative velocity v between projectile and target is assumed to be large. 

The first step in the eikonal approximation is to assume the second-order deriva- 
tive P^/ljx negligible with respect to the first-order derivative vPz, 

— f(^,R)«vPzf(<^,R). (4.3) 
This first step leads to the second-order equation (but now first-order in Z), 

ihv^W{^,h,Z) = [//o-£o + Vpr(^,R)]'?(^,b,Z), (4.4) 

where the dependence of the wave function on the longitudinal Z and transverse 
b parts of the projectile-target coordinate R has been made explicit. This equation 
is mathematically equivalent to the time-dependent Schrodinger equation ( 13.1b for 
straight-line trajectories with t replaced by Z/v. It can thus be solved using any of 
the algorithms cited in Sec. 13.41 However, contrary to time-dependent models, it is 
obtained without the semiclassical approximation. The projectile-target coordinate 
components b and Z are thus quantal variables. Interferences between solutions ob- 
tained at different b values are thus taken here into account. This first step is known 
as the dynamical eikonal approximation (DEA) lfT2l [T3l . 
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The transition matrix element for elastic scattering into direction Q = {9 ,(p) of the 
final momentum K = {K,Q) of the projectile in the cm. frame reads ll37ll 



Tfi = (e''^-'*0':''°(£o,^)|Vpr(^,R)|f (^'''(^,R)), (4.5) 

where ^''^'^o) t^e exact solution of the Schrodinger equation ( I2.4l l with the asymp- 
totic condition (12.6b . By using Eqs. (14. 11 1. ( 12.1b . and ( 14.41 ). one obtains the approxi- 
mation lfT2l 

« ihv I ^/Re-'i''^(</.^r"(^o,^)|'f'(^«)(^,R)), (4.6) 

where the transfered momentum q = K — KZ is assumed to be purely transverse, 
i.e. exp[/(K • Z — K)], is neglected. The norm of q is linked to the scattering angle by 

q = 2Ksme/2. (4.7) 

Let us define the elastic amplitude 

^Si(b) = (0r"(£o,^)|^'f^''n^,b,Z^ +^)) - 5^,^^. (4.8) 
The transition matrix element ( 14.6b reads after integration over Z, 

Tfi = ihv [ dbe-''i-^e'^^o-M'o)<PRs^J^^, (bX), (4.9) 



where (pu is the azimuthal angle characterizing b. The phase factor exp[/(Mo — 
Mq)(Pr] arises from the rotation of the wave functions when the orientation of b 
varies (131 . The integral over (pR can be performed analytically, which leads to the 
following expression for the elastic differential cross section llT3l 



(4.10) 



where J„, is a Bessel function 1*291. From Eq. ( 14. 10b . one can see that contrary to the 
semiclassical approximation ( 13. 5b . the eikonal elastic cross section is obtained as a 
coherent sum of elastic amplitudes over all b values. This illustrates that quantum 
interferences are taken into account in the eikonal framework. 
The transition matrix element for dissociation reads 

Tf, - (e''^'«0[-'(£,<^)|ypH<^,R)|f '''"^(<^,R)), (4.11) 
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where K' = {K' ,Q) is the final projectile-target wave vector. One can then proceed 
as for the elastic scattering. Using Eqs. ( 14.1b . i2.2i . and ( 14.41 ). taking into account 
the energy conservation. 



2 1^2 



2^ 



2vl2 



-Eo = 



2jx 



and assuming the transfered momentum q 
transition matrix element is expressed as 



K' 



£, (4.12) 
KTj to be purely transverse, the 



Tfi « ihv J dhe-"i-^sl''°^ 



-"bu 



with the breakup amplitude 



'(£,^)|'f'(^«)(^,b,Z- 



The differential cross section for breakup is given by 



da 



dk^dEdQ 2Jq 



— E 



dhe 



-/q-br.('Wo) 
'^bu 



(4.13) 



(4.14) 



(4.15) 



where the proportionality factor depends on the phase space. Like the elastic scat- 
tering cross section ( 14.10b . it is obtained from a coherent sum of breakup amplitudes 
(14.14b . confirming the quantum-mechanical character of the eikonal approximation. 
Here also, the integral over (j?^ can be performed analytically and leads to Bessel 
functions lfT3l . 

By integrating expression ( 14. 15b over unmeasured quantities, one can obtain the 
breakup cross sections with respect to the desired variables, like the internal excita- 
tion energy of the projectile. Since these operations depend on the projectile internal 
structure, we delay the presentation of some detailed expressions to Sees. |6] and |7] 
treating of two-body lfT3l and three-body ||T9| breakup. 



4.3 Standard eikonal approximation 

In most references, the concept of eikonal approximation involves a further simpli- 
fication to the DEA 1381 [30| . This adiabatic, or sudden, approximation consists in 
neglecting the excitation energy of the projectile compared to the incident kinetic 
energy. It comes down to assume the low-lying spectrum of the projectile to be de- 
generate with its ground state, i.e. to consider the internal coordinates of the projec- 
tile as frozen during the reaction |30|. This approximation therefore holds only for 
high-energy collisions that occur during a very brief time. This second assumption 
leads to neglect the term Hq — Eq in the DEA equation (14.41 ) which then reads 
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ihv-^^i^M ^VpTi^,R)'Pi^,h,Z). (4.16) 

The solution of Eq. ( I4.16l l that follows the asymptotic condition ( 12.61 ) exhibits the 
well-known eikonal form lfm[36l 



$^(^'')(^,b,Z) =exp 



<«'"»(£o,^). (4.17) 



After the collision, the whole information about the change in the projectile wave 
function is thus contained in the phase shift 

Z(s^,b) = - ^ £yPT{^,R)dZ. (4.18) 

Due to translation invariance, this eikonal phase x depends only on the transverse 
components b of the projectile-target coordinate R and of the projectile inter- 
nal coordinates Cross sections within this standard eikonal approximation are 
obtained as explained in Sec. 14.21 replacing 

Being obtained from the adiabatic approximation, expressions ( 14.17b and (14.18b 
are valid only for short-range potentials. For the Coulomb interaction, the assump- 
tion that the reaction takes place in a short time no longer holds, due to its infi- 
nite range. The adiabatic approximation thus fails for Coulomb-dominated reactions 
ll30l . Besides imprecise uses of a cutoff at large impact parameters ||39l , there are 
two ways to avoid this problem. The first is not to make the adiabatic approxima- 
tion, i.e. to resort to the more complicated DEA (see Sec. 14.1b . The second is to 
correct the eikonal phase for the Coulomb interaction as suggested in Ref. |[T4| (see 
Sec. |4.4b . Nevertheless, as shown in Ref. fTJI , the Coulomb divergence does not 
affect eikonal calculations performed on light targets at high enough energies. Most 
of the nuclear-dominated reactions can thus be analyzed within an eikonal model 
including the adiabatic approximation (see, e.g., Ref. 00). 

Fig.[2]illustrates the difference between the DEA (full line), the usual eikonal ap- 
proximation (dashed line) and the semiclassical approximation (dotted line) when 
Coulomb dominates. It shows the breakup cross section of ''Be on Pb at 69 
MeV/nucleon for a '"Be-n relative energy of 0.3 MeV as a function of the P-T 
scattering angle. As explained above, the usual eikonal approximation diverges for 
the Coulomb-dominated breakup, i.e. at forward angles. The DEA, which does not 
include the adiabatic approximation, exhibits a regular behaviour at these angles. 
Interestingly, the semiclassical approximation follows the general behaviour of the 
DEA, except for the oscillations due to quantum interferences between different 
b values. The DEA has therefore the advantage of being valid for describing any 
breakup observable on both light and heavy targets. 

The nuclei studied through breakup reactions being exotic, it may be difficult, if 
not impossible, to find optical potentials that describe the scattering of the clusters 
by the target. One way to circumvent that problem is to resort to what is usually 
known as the Glauber model B2l [38l [30l [36l . This model has been mostly used 
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Fig. 2 Breakup cross section of " Be on Pb at 69 MeV/nucleon as a function of the P-T scattering 
angle in the P-T cm. frame for a '"Be-n energy E — 0.3 MeV. Calculations are performed within 
the DEA, usual eikonal, and semiclassical approximations 1411 . 

to calculate total and reaction cross sections. At the optical-limit approximation 
(OLA) of the Glauber model, correlations in the cluster and target wave functions 
are neglected. The nuclear component of the eikonal phase shift for cluster / is then 
expressed as a function of the densities pr of the target and Pi of the cluster, and of 
a profile function 1 — c'^nn that corresponds to an effective nucleon-nucleon inter- 
action. The nuclear component of the eikonal phase shift is approximated by ll30l 



where St and s, are the transverse components of the internal coordinates rj- of the 
target and r, of cluster /, respectively, and b, is the transverse component of the 
cm. coordinate of cluster /. The OLA is therefore equivalent to the double-folding 
of an effective nucleon-nucleon interaction. The density of the target can usually 
be obtained from experimental data. The cluster density being unknown, it has to 
be estimated from some structure model, like a mean-field calculation. The profile 
function is usually parametrized as 1301 l43l 



where is the total cross section for the N-N collision, aNN is the ratio of the 
real part to the imaginary part of the N-N scattering amphtude, and J3nn is the slope 
parameter of the N-N elastic differential cross section. These parameters depend 
on the nucleon type (p or n) and on the incident energy. Their values can be found 
in the literature (see, e.g., Ref. ||43l ). The validity of the Glauber approximation is 
discussed in Ref. 1441 . 




(4.19) 




(4.20) 
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4.4 Coulomb-corrected eikonal approximation 

The eikonal approximation gives excellent results for nuclear-dominated reactions 
ll30l [T3I . However, as mentioned above, it suffers from a divergence problem when 
the Coulomb interaction becomes significant. To explain this, let us divide the 
eikonal phase ( 14.181 ) into its Coulomb and nuclear contributions 

Z(s^,b)=;t^^(/7) + zC(s^,b) + z^(s^,b). (4.21) 

In this expression, Xpj '^^e global elastic Coulomb eikonal phase between the 
projectile and the target. However, Coulomb forces not only act globally on the 
projectile, they also induce 'tidal' effects due to their different actions on the var- 
ious clusters. The tidal Coulomb phase is due to the difference between the 
cluster-target and projectile-target bare Coulomb interactions. The remaining phase 
x'^ contains effects of the nuclear forces as well as of differences between Coulomb 
forces taking the finite size of the clusters into account and the bare Coulomb forces. 

At the eikonal approximation, the integral ( |4.18t defining Xpt diverges and must 
be calculated with a cutoff ifTTllSOl . Up to an additional cutoff-dependent term that 
plays no role in the cross sections, it can be written as ll36l 

X^PT{b)^2iq\n{Kb), (4.22) 

where appears the projectile-target Sommerfeld parameter Tj defined in Eq. ( 12.8b . 
The phase ( I4.22l i depends only on b. 

The tidal Coulomb phase is computed with Eq. ( 14.181 1 for the difference between 
the bare Coulomb interactions for the clusters in the projectile and the global P-T 
Coulomb interaction. 

It can be expressed analytically. Because of the long range of the El component 
of the Coulomb force, this phase behaves as l/b at large distances lfT3l[T6l . In the 
calculation of the breakup cross sections ( 14.15b . the integration over bdb diverges for 
small q values, i.e. at forward angles, because of the corresponding 1 /b asymptotic 
behaviour of the breakup amplitude, as illustrated in Fig.|2] This divergence occurs 
only in the first-order term ix'' of the expansion of the eikonal Coulomb amplitude 
exp(;x'^). 

As seen in Sec. 13.31 the first order approximation ( 13.151 ) decreases exponentially 
at large b and hence does not display such a divergence. A plausible correction is 
therefore to replace the exponential of the eikonal phase according to lfT4l[T5l 

where x^^ is the result of first-order perturbation theory ( 13.10b . 



Breakup reaction models for two- and three-cluster projectiles 



17 



Note that because of the phase e""^/*', the integrand in Eq. ( 14.251 ) does not exhibit 
a translational invariance. The first-order phase x^'^ depends on all internal coordi- 
nates of the projectile. When the adiabatic approximation is applied to Eq. ( 14.251 ), 
i.e. when (O is set to 0, one recovers exactly the Coulomb eikonal phase ( 14.231 ). 
This suggests that without adiabatic approximation the first-order term of exp{ix^) 
would be /j^^ ( |4.25l l instead of i^'^ ( |4.23t , intuitively validating the correction 
( 14.241 ). Furthermore, since a simple analytic expression is available for each of the 
Coulomb multipoles (see Sec. 13.31 ). this correction is easy to implement. 

With this Coulomb correction, the breakup of loosely-bound projectiles can be 
described within the eikonal approximation taking on (nearly) the same footing both 
Coulomb and nuclear interactions at all orders. This approximation has been tested 
and validated for a two-body projectile in Ref. lT6l . Note that in all practical cases 
lfT5l[T6l[T9l . only the dipole term of the first-order expansion ( 13.131 ) is retained to 
evaluate x'''^ 

Fig. [3] illustrates the accuracy of the CCE for the breakup of "Be on lead at 
69 MeV/nucleon lfT6l . The figure presents the parallel-momentum distribution be- 
tween the '''Be core and the halo neutron after dissociation. This observable has 
been computed within the DEA (full line), which serves as a reference calculation, 
the CCE (dotted line), the eikonal approximation including the adiabatic approx- 
imation (dashed line), and the first-order of the perturbation theory (dash-dotted 
line). The usual eikonal approximation requires a cutoff at large impact parameter 
to avoid divergence. The value bmax =71 fm is chosen from the value prescribed 
in Ref. ||39l . At the first order or the perturbation theory, the nuclear interaction is 
simulated by an impact parameter cutoff at b^in = 15 fm. 
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Fig. 3 Breakup of "Be on Pb at 69 MeV/nucleon. The parallel-moinentuiTi distribution between 
the "'Be core and the halo neutron is computed within the DEA, the CCE, the eikonal approxi- 
mation including the adiabatic approxiination, and the first-order of the perturbation theory [16]. 
Reprinted figure with perinission from Ref. [HI. Copyright (2008) by the Amer ican Physical So- 
ciety. 
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We first see that the magnitude of the CCE cross section is close to the DEA one, 
whereas, the other two approximations give too large (eikonal) or too small (first 
order) cross sections. Moreover the CCE reproduces nearly perfectly the shape of 
the DEA distribution. In particular the asymmetry, due to dynamical effects, is well 
reproduced. This result suggests that in addition to solving the Coulomb divergence 
problem introduced by the adiabatic approximation, the CCE also restores some 
dynamical and higher-order effects missing in its ingredients, the usual eikonal ap- 
proximation and the first order of the perturbation theory. 



5 Continuum-discretized coupled-channel method 

The CDCC method is a fully quantal approximation which does not imply some 
restriction on energies. Its main interest lies in low energies where the previous 
methods are not valid. The principle of the CDCC method is to determine, as ac- 
curately as possible, the scattering and dissociation cross sections of a nucleus with 
a simplified treatment of the final projectile continuum states. To this end, these 
states describing the relative motions of the unbound fragments are approximately 
described by square-integrable wave functions at discrete energies. The relative mo- 
tion between the projectile and target and various cross sections can then be obtained 
by solving a system of coupled-channel equations. The number of these equations 
and hence the difficulty of the numerical treatment increase with increasing energy. 

The CDCC method was suggested by Rawitscher j45| and first applied to 
deuteron + nucleus elastic scattering and breakup reactions. It was then extensively 
developed and used by several groups Inl |46l ll?] |48l |49] |50l |5l]. Its interest has 
been revived by the availabihty of radioactive beams of weakly bound nuclei disso- 
ciating into two |[l8l|47l|48l|451|50l|5ll or three fSS] |53l |2 2 1 1 fragments. 

We assume that the breakup process leads to clusters and that the cluster- 
target interactions do not depend on the target spin. The projectile wave func- 
tions describing A^-body bound states at negative energies E'^g and 
(^^*'(£',^) describing A^-body scattering states at positive energies E are defined 
with Eq. ( 12. 11 1. Since the total angular momentum of the projectile-target system 
is a good quantum number, the first step consists in determining partial waves of 
the (A^+ l)-body Hamiltonian ( 12.4b . The general partial wave function for a total 
angular momentum Jj can be expanded over the projectile eigenstates as 



In this expansion, index B runs over the bound states of the projectile. The total 
angular momentum Jj results from the coupling of the orbital momentum L of the 
projectile-target relative motion with the total angular momentum / of the projectile 



ZJt B 




(5.1) 
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State. The relative-motion partial waves ij/j^-g and i/y^^ are unknown and must be 
determined. The parity is given by the product of (—1)^ and the parity of . 
The first term of Eq. ( 15.1b represents the elastic and inelastic channels while the 
second term is associated with the breakup contribution. However, the presence of 
the continuum renders this expression intractable. 

The basic idea of the CDCC method is to replace wave function ( 15.1b by 

^7xMt(r^^) = ^ v/j-,„(R)]^TMx^ (5.2) 

LJzn 

where the functions ^(ff{^) = ^^*'(£'4i^) represent either bound states (E-^g < 
0) or square-integrable approximations of continuum wave functions {E'^^ > 0) at 
discrete energies 

Ei„ = {^i^{^)m^i^{^)). (5.3) 

Approximation ( 15.21 ) resembles usual coupled-channel expansions and can be treated 
in a similar way. 

In practice, two methods are available to perform the continuum discretization. In 
the "pseudostate" approach, the Schrodinger equation ( 12.1b is solved approximately 
by diagonalizing the projectile Hamiltonian Hq either within a finite basis of square- 
integrable functions or in a finite region of space. In both cases, square-integrable 
pseudostates are obtained. This approach is simple but there is little control on 
the obtained energies E'^^. Therefore, it is customary to keep only the pseudostates 
with energies below some limit Em^x- 

The alternative is to separate the integral over E in Eq. ( 15.11 ) into a limited number 
of small intervals, or "bins", [En-i,E„] which may depend on J and to use in each 
of them some average of the exact scattering states in this range of energies |[T7ll45l 
|46l|42l. This "bin" method provides the square-integrable basis functions 

= 7^ r <^f{E.i)UE)dE, (5.4) 

where the weight functions /„ may also depend on J. Such states are orthogonal 
because of the orthogonaUty of the scattering states and they are normed if W,, is the 
norm of /„ over \,E„]. Using Eq. ( 15.4b . their energy ( 15.31 ) is given by 

< = i C" \fn{E)\'EdE. (5.5) 

Here also, a maximum energy fimax = E„^^,^^ is chosen. In practice, these basis states 
are usually constructed by averaging the scattering states 0^^^ normahzed over 
the wave number k, often within equal momentum intervals ifTSl . 
The total wave function (15.21 ) can be rewritten as 



(5.6) 
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where c represents the channel LJxn and 

By inserting expansion (15.6b in the Schrodinger equation (I2.4l i and using Eq. ( 15.3b . 
the relative wave functions up^ are given by a set of coupled equations 

«?W+L^/M^)«?W==0, (5.8) 

where = E(„. The sum over L is truncated at some value Lmax. The sum over 
the pseudo-states or bins is limited by the selected maximum energy fimax. The 
CDCC problem is therefore equivalent to a system of coupled equations where the 
potentials are given by 

<c'(^) - {i>i'''HnR,^)\VPT{R,^)\^^J''^nR,^)). (5.9) 

This matrix element involves a multidimensional integral over and over the in- 
ternal coordinates ^ . In general, the potentials are expanded into multipoles corre- 
sponding to the total angular momentum operator Jj of the system. This may allow 
an analytical treatment of angular integrals. 

System (15.8b must be solved with the boundary condition 

u'j{R) —> v-'/2 [l,{K,R)5,c, - Oc{KcR)S'jJ , (5.10) 

R — ^co 

where cq is the incoming channel. The asymptotic momentum in channel c reads 

K, = ^2pL{Ej-E,)/n^, (5.11) 

and \'c = TiKc/ jJ. is the corresponding velocity. In Eq. ( 15.10b . Ic = Gc — iFc and Oc ~ 
I* are the incoming and outgoing Coulomb functions, respectively ll29l . and the 
element ^cJ,, of the collision matrix is the amplitude for populating channel c from 
initial channel cq. 

Various methods have been developed to solve system ( 15.81 ) (see, e.g., Ref. 1541 ). 
A convenient approach is the /^-matrix formalism ll55l . which is both simple and 
accurate. The configuration space is divided into two regions: the internal {R < a) 
and external (R > a) regions, where a is the channel radius. In the external region, 
the potential matrix defined by Eq. ( 15. 9b can be well approximated by its diagonal 
Coulomb asymptotic form. Hence the wave function is replaced by combinations 
of Coulomb functions. In the internal region, the radial wave functions u'jy can be 
expanded over some basis 1551 . A significant simplification occurs when using La- 
grange functions l56l l57l ISTl . 

A scattering wave function verifying the initial condition ( 12.6b is then constructed 
with the different partial waves. Inserting this CDCC approximate wave function in 
Eq. ( 14.111 ) enables calculating transition matrix elements towards pseudostates or 
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Fig. 4 'Be angular distribution after the breakup of **B on '^'Ni at 25.8 MeV computed within a 
CDCC model 1181 . The convergence of the numerical scheme is illustrated with vaiious maximum 
values of the ^Be-p relative orbital momentum / in the continuum and various maximum values of 
multipole order A of the potential expansion in Eq. )5.9t : / < 3, A < 2 (dashed line), / < 3, A < 3 
(full line), / < 4, A < 4 (dash-dotted line). Experimental data from Ref. |62'|. Reprinted figure with 
permission from Ref. 1181 . Copyright (2001) by the American Physical Society. 

bin states as a function of the collision matrices 5'''t (see Eq. (5) of Ref. IfTSlI ). Since 
these transition matrix elements are obtained only at discrete energies E^^, they must 
be interpolated in order to obtain breakup cross sections at all energies. 

The CDCC method has first been applied to two-body projectiles. As an example, 
Fig. |4]shows the convergence of the breakup of **B on ^^Ni at 25.8 MeV. The con- 
vergence concerns the set of partial waves / in the ^Be-p continuum of the projectile 
and the number of multipoles in the expansion of the potential appearing in matrix 
elements ( 15.9b . The validity of CDCC has been tested for breakup observables in 
a comparison with three-body Faddeev calculations ||24ll • The agreement between 
both sets of results is good except when the coupling with the transfer channel is 
important. 

Let us also mention extensions beyond the simple two-body model of the projec- 
tile by allowing the core to be in an excited state 15811591 . These references present 
total cross sections for the breakup on a ^Be target of "Be into "^Be + n and of '^C 
into '^c + n calculated by including core deformations. This extension of CDCC 
known as XCDCC leads to very long computational times. 

The extension of CDCC to three-body projectiles is more recent ll52l l60l l53l 
I20II2T]. The calculations are still much more time-consuming since the projectile 
wave functions are much more complicated (see Sec. |7]i. Consequently, the cal- 
culation of the potential matrix elements i5.9i raises important numerical difficul- 
ties. At present, converged calculations are mainly restricted to elastic scattering 
II52I l53l l60l . Most breakup calculations still involve limited bases and/or simplify- 
ing assumptions 1201 12TI but these Umitations can be overcome [61 1. 




22 



D. Baye and P. Capel 



6 Breakup reactions of two-body projectiles 
6.1 Two-clustermodel 

Most of the reaction models have been applied assuming a two-cluster structure of 
the projectile. In this section, we specify the expression of the internal Hamiltonian 
of the projectile and the set of coordinates usually considered in practical appli- 
cations. We then illustrate the models presented in Sec. |2] and the approximations 
explained in Sees. [3H5] with various applications to the study of exotic nuclei and 
nuclear astrophysics. 

We consider here projectiles made up of a single fragment / of mass nif and 
charge Zfe, initially bound to a core c of mass nic and charge Z^e. The core and 
fragment are assumed to have spins and Sf. The internal structure of these clus- 
ters and of the target is usually neglected although some structure effects can be 
simulated by the effective potentials. 

Let us now particularize the general formalism (I2.1b -( l23] l to the present case. 
The internal coordinates ^ represent the relative coordinate r = r/^ — . The struc- 
ture of the projectile is described by the two-body internal Hamiltonian 

2-llcf 

where /i^/ = nicinf/mp is the reduced mass of the core-fragment pair (with mp = 
nic+mf), p is the momentum operator of the relative motion and Vcf is the potential 
describing the core-fragment interaction. This potential usually includes a central 
part and a spin-orbit coupling term in addition to a Coulomb potential. In many 
cases, the potential is deep enough to contain unphysical bound states below the 
ground state. These unphysical or forbidden states are useful because they allow 
the wave function representing the physical ground state to exhibit the number of 
nodes expected from the Pauli principle, as obtained in microscopic descriptions 
ll63l . Although these forbidden states do not play any role in the core-fragment 
scattering, they could affect breakup properties. However, as shown in Ref. [^3, 
their presence can be ignored because their effect is negligible. 

Let k be the wave vector describing the asymptotic relative motion between 
the fragments in the projectile continuum. The corresponding energy is thus E = 
Tp'k^/ljicf- Notation t in Eq. (12.1b corresponds here to the coupling mode, i.e. to 
the total spin S of the projectile and the relative orbital momentum /. The wave 
functions defined in Eq. ( 12.11 ) read 

(£ , r) = r- ' /' [Yi {Q ) ® XsV^ufs {k, r) , (6.2) 

where Xs is a spinor resulting from the coupling of Sc and Sf. The radial functions 
u'ig{k,r) are normalized according to {u'l^{k, r)\u'jg{k' , r)) = 5{k — k'). The notation 
in Eq. ( 12.2b represents the direction £1]^ of k and the spin orientations Vc and V/ of 
the core and fragment spins Sc and Sf. Relation (12.3b between continuum eigenstates 
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of Hq becomes 

Cv.v,(£.r) = i ^ {s,SfV,Vf\Sv){lSM-vv\JM)Yf'-'*{£l,)<^fs''{E,r) (6.3) 

ISJM 

with the property ((/)|2;,'^^^^,(£,r)|(/)<27'^,^, (£',r)) = 5(k-k')Vv^5v^v}- Notice that 

notations T and are model dependent and would be quite different if a tensor 
interaction were included in Vcf- A detailed description of the simple case Sc—Sf — 
can be found in Ref . Il26l . 




Fig. 5 Jacobi set of coordinates: 

r is the projectile internal coordinate, and 

R = b + ZZ is the target-projectile coordinate. 



Within this framework the description of the reaction reduces to the resolution of 
a three-body Schrodinger equation (12.41 ) that reads, in the Jacobi set of coordinates 
illustrated in Fig.|5] 



yp2 

— +//o + y/.r(r,R) 



f(r,R) =£Tf(r,R). 
The projectile-target interaction (12. 5t then reads 



VfT R+— r 
* mp 



ypr(R,r) = y,r R-^r 

mp 



(6.4) 



(6.5) 



where Vct and Vfp are optical potentials that simulate the core-target and fragment- 
target interactions, respectively. 

For a two-body projectile, the DBA breakup cross section ( 14.15b becomes Eq. (46) 
of Ref. IIT3I . Integration over Q/^ and summation over Vc and Vf lead to the energy 
and angular distribution of the fragments in the P-T cm. rest frame. With the nor- 
malization of the positive-energy states given above, it reads lfT3l 



dEdQ 



KK' \ - \ - 



bdbJ\M-M,Mb)S^umib) 



(6.6) 
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where S^if^jj are coefficients of a partial-wave expansion of the breakup ampUtude 
(14.141 1 (see Eq. (43) of Ref. lfT3l ). Breakup cross sections are mainly expressed as 
energy distributions da\m/dE as a function of the energy of the relative motion be- 
tween the fragments. They are obtained by integrating ( 16.6b over Q. However, most 
experimental data concern angular distributions or distributions of the core momen- 
tum in the laboratory frame. Note that, in addition, theoretical results should be 
convoluted with the experimental acceptance and resolution. A change of frame for 
the theoretical results is thus in general not sufficient to allow a fruitful comparison 
with experiment. 



6.2 Two-body breakup of exotic nuclei 

A first information that one can extract from experiment concerns the separation 
energy of the halo neutrons. Indeed, the shape of the breakup cross section and, in 
particular, its maximum are sensitive to this energy as can be shown at first order 
of perturbation theory with rather simple models based on the asymptotic behaviour 
of the halo wave function ||65]| . An example is given by the breakup of ^^C on lead 
at 67 MeV/nucleon. In Fig.|6] a semi-classical calculation with a ^^C + n two-body 
model shows that the shape of the experimental data is much better reproduced if 
the binding energy of '^C is raised from the recommended value 0.53 MeV to 0.65 
MeV Il34l, 




"C + Pb -> '°C + n + ™Pb 
Ek„_ = 67 MeV/nucleon 

Coul. dynamical 

Coul. + Nucl. dynamical 

Nucl. dynamical 
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Fig. 6 Breakup of on Pb at 67 IVIeV/nucleon: semi-classical cross sections for two different 
binding energies of the projectile: 0.53 MeV (upper panel) and 0.65 MeV (lower panel) 1341 . 
Experimental data from Ref. 1661 . Reprinted figure with permission from Ref. (34). Copyright 
(2001) by the American Physical Society. 
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Fig. 7 Breakup of " Be on a C target at 67 MeV/nucleon: calculation performed in a semi-classical 
model 1691 . Experimental data from Ref. 1351 . 

Indirect information can also be obtained on the spin of the ground state of the 
halo nucleus when few rather different orbital momenta are probable. The magni- 
tude of the cross section is very sensitive to the orbital momentum / of the ground 
state. A study of the one-neutron removal cross section from ^^Ne described in a 
simple ^"Ne + n model allows to rule out the prediction 7/2^ of the naive shell 
model and to confirm the value 3/2^ resulting from a shell inversion Ii67ii68j . 

Nuclear-induced two-body breakup on light targets is an interesting tool to 
observe resonances of a halo nucleus and to assess some of their properties. In 
Fig. I?] are displayed experimental data on the ''Be breakup on a C target at 
67 MeV/nucleon ||35l . These data present a broad bump near the location of a known 
resonance with an assumed spin-parity 5/2+. The bump width is however broader 
than the known resonance width. A semi-classical calculation (dashed line) based 
on a '"Be + n model reproduces the shape of the data very well after convolution 
with the experimental resolution (full line). Moreover the d5/2 component of the 
theoretical cross section (dotted hne) resonates and confirms the 5/2+ attribution. 

Breakup reactions are also used to infer the spectroscopic factor of the dominant 
configuration in the core-nnucleon structure of halo nuclei 1251 l40l . Various theo- 
retical studies have been performed to assess the sensitivity of breakup calculations 
to the projectile description ITOIITTI . These studies have revealed that the breakup 
cross sections not only depend on the initial bound state of the projectile, but are 
also sensitive to the description of its continuum [70J . Moreover it has been shown 
that, for loosely-bound projectiles, only the tail of the wave function is probed in 
the breakup process and not its whole range 17111721 . These studies therefore sug- 
gest that one should proceed with caution when extracting spectroscopic factors of 
weakly -bound nuclei from breakup measurements, as other structure properties, like 
the continuum description, may hinder that extraction. 

As mentioned earlier, many Coulomb-breakup experiments have been analyzed 
within the framework of the first order of the perturbation theory (see Sec. 13.3b . 
In order to assess the validity of that approximation, various authors have com- 
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Fig. 8 Influence of the couplings inside the continuum 1 28) . Time evolution of the numerical 
breakup probability per energy unit l |3.7) for "Be impinging on Pb at about 45 MeV/nucleon 
for a '"Be-n relative energy £ = 1.5 MeV and an impact parameter b = 100 fm. Reprinted figure 
with permission from Ref. 1281 . Copyright (2005) by the American Physical Society. 



pared perturbation calculations to numerical resolutions of the time-dependent 
Schrodinger equation 17311741 [75] |28l . These studies have shown that, in many cases, 
breakup cannot be modelled as a one-step process from the initial bound state to- 
wards the continuum and that higher-order effects should be considered for a reliable 
description of the reaction. In particular, they indicate that significant couplings are 
at play inside the continuum. To illustrate this, Fig.|8]displays the time evolution of 
the breakup probability per energy unit ( 13.7b for the collision of ^'Be on Pb at about 
45 MeV/nucleon computed within the time-dependent '"Be + n model of Ref. lITOl . 
The obtained value is divided by its evaluation at the first-order of the perturba- 
tion theory ( 13.15b at f ^ +°o. After a sharp increase at the time of closest approach 
t = 0, the breakup probability (full line) oscillates and then stabilizes at a value 
which differs by about 5% from its first-order estimate. Although the total breakup 
probability becomes stable, its partial-wave composition still varies: the dominant 
p wave contribution (dash-dotted line) is depleted towards the s (dotted line) and 
especially d (dashed line) ones. This signals couplings inside the continuum, which 
may affect the evaluation of breakup observables 1281 l75l . We will see in the next 
section that it may perturb the analysis of breakup reactions of astrophysical interest 
Il76i|77l|78l|79l. 



6.3 Application to nuclear astrophysics 

Radiative-capture reactions are a crucial ingredient in the determination of the reac- 
tion rates in nuclear astrophysics. However the difficulty of their measurement and, 
in some cases, the scatter of the results has raised interest in indirect methods where 
the time-reversed reaction is simulated by virtual photons in the Coulomb field of 
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a heavy nucleus ISOllSTl . The radiative-capture cross section can be extracted from 
breakup cross sections if one assumes that the dissociation is due to El virtual pho- 
tons and occurs in a single step. A typical example is the ^Be(p,7)^B reaction which 
has been studied with the breakup of into ^Be-np on heavy targets at different 
energies ||82ll83l|84l[85l|62l[H6ll821l88l. 

Though appealing, the breakup method also faces a number of difficulties. First, 
while many reactions are dominated by an El transition, an E2 contribution to the 
breakup cross section may not be negligible 1731 . Second, higher-order effects, i.e. 
transitions from the initial bound state into the continuum through several steps 
may not be negligible 1731 l74l l75l |28l . Finally, the nuclear interactions between 
the projectile and the target may interfere with the Coulomb interaction l34l [TOll . 
Therefore elaborate reaction theories must be used to interpret the experimental 
data. 
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Fig. 9 ^'B Coulomb breakup on Pb at 44 MeV/nucleon. Parallel-momentum distribution of the 
Be core corresponding to various angular cuts calculated in a DBA model (TT)- (a) Influence of 
nuclear and Coulomb interactions on the calculation, (b) Effects of the various multipoles of the 
Coulomb interaction, (c) Role of the higher-order effects on the calculation. Experimental data 
from Ref. (84] ■ 

The experiments on the breakup of ^B have been analyzed in a number of papers 
lfT8ll89ll90ll9ni76ll77l . Figure|9]shows a comparison between the experimental data 
of Ref. ll84l and DEA calculations fTT l. Without adjustable parameters, the calcula- 
tions (full lines) fairly reproduce the asymmetry exhibited by the data which could 
not be well explained in earlier works l89l |9T1 . The three panels of Fig. |9]illustrate 
the influence of various approximations upon the calculation fTTll. The left panel 
illustrates that nuclear P-T interactions can be neglected when data are restricted to 
forward angles. The central panel confronts a dynamical calculation including only 
the dipole term of the Coulomb interaction (dashed lines) to the full calculation, 
indicating that higher multipoles have a significant effect on the breakup process. 
The right panel compares the dynamical calculation to its first-order approxima- 
tion (dot-dashed lines), emphasizing the necessity to include higher-order effects 
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Fig. 10 Breakup of '^C on Pb at 68 MeV/nucleon. The experimental energy distribution measured 
for two scattering-angle cuts |98| is confronted to the time-dependent calculation of Ref. |79| . 
Reprinted figure with permission from Ref. 1791 . Copyright (2009) by the American Physical So- 
ciety. 



in breakup calculations. These results show that some of the assumptions of the 
breakup method ISOl [STj are not valid. It is therefore difficult to infer the accuracy 
of the S factors extracted from breakup cross sections. 

An interesting problem was raised by the ''^C(n,7)'^C capture reaction. The mea- 
sured cross sections for the Coulomb breakup of '^C 19211931 provided an S factor 
which disagreed with direct measurements ll94ll95]| . Moreover, theoretical analyses 
indicated that the Coulomb-breakup cross sections were inconsistent with informa- 
tion obtained from '^F by charge symmetry and with microscopic models ||96l . A 
new measurement 1971 |98l has obtained breakup cross sections that fully agree with 
properties of the mirror system and with theory 17811791 . These theoretical analyses 
show that a fully dynamical calculation, taking proper account of higher-order ef- 
fects is necessary to correctly analyze the breakup measurements, in agreement with 
the analysis of the Coulomb breakup of Ref. f771 . They also indicate that includ- 
ing both Coulomb and nuclear interactions as well as their interferences is neces- 
sary to correctly reproduce data at large scattering angles. In this way a very good 
agreement can be obtained between direct and indirect measurements of the S factor. 
Fig.[TO]displays the breakup cross section of '^C on Pb measured at 68 MeV/nucleon 
l98il and its comparison to the theoretical calculation of the time-dependent model 
of Ref. l79l . The dotted lines show the direct results of the calculation, while the full 
lines correspond to these results folded by the experimental resolution and scaled to 
the data. 



Breakup reaction models for two- and three-cluster projectiles 



29 



7 Breakup reactions of three-body projectiles 
7.1 Three-cluster model of projectile 

Let us consider a system of three particles, the core with coordinate r^, mass nic and 
charge Z^e and two fragments with coordinates ri and r2, masses m\ and m2, and 
charges Z\e and Z2e. The projectile mass is nip = iric + mi2 with mi2 = mi +m2. 
After removal of the cm. kinetic energy Tc.m., the three-body Hamiltonian of this 
system in Eq. ( 12. Il l can be written as 



where Vij is an effective potential between particles / and j ~ c, 1,2). We as- 
sume that these interactions involve central, spin-orbit and Coulomb terms. These 
potentials may contain unphysical bound states below the two-cluster ground state 
to simulate effects of the Pauli principle. These forbidden states must be eliminated 
from the three-body wave functions either with pseudopotentials ll99l or with super- 
symmetric transformations HI 001 IIOII . 

Various resolution techniques can be considered for obtaining the wave functions 
of a three-body projectile. A first option is to describe this projectile with an expan- 
sion in Gaussian functions depending on Jacobi coordinates Ill02ll30ll20l . For bound 
states, the wave functions can be obtained from a variational calculation. Well es- 
tablished techniques allow systematic calculations of the matrix elements Ill02|[30l . 
Calculations are then simpler when the interactions are expressed in terms of Gaus- 
sians. At negative energies, this type of expansion may however have convergence 
problems in the description of extended halos. At positive energies, it is convenient 
to obtain pseudostates but not convenient to obtain scattering states. 

Let us describe another efficient tool to deal with three-body systems, the formal- 
ism of hyperspherical coordinates. It is especially interesting when the two-cluster 
subsystems are unbound so that only a three-body continuum exists. Notation ^ of 
Sec.|2]represents here five angular variables and one coordinate with the dimension 
of a length, the hyperradius (see Refs. 11031 11041 1105II for details). Four angular 
variables correspond to physical angles and the fifth one is related to a ratio of co- 
ordinates defined below in Eq. ( 17.6b . The wave functions are expanded in series of 
hyperspherical harmonics, i.e. a well known complete set of orthonormal functions 
of the five angular variables. The coefficients are functions of the hyperradius and 
can be obtained from variational calculations. Scattering states can be obtained from 
extensions of the R matrix theory 11061 [55]| . A drawback of this method is that the 
hyperspherical expansion may converge rather slowly. 

With the dimensionless reduced masses jJ-c^n) = mcm\2/mpmN and pLi2 = mim2 
/m\2mN where m^ is the nucleon mass for example, the internal coordinates ^ are 
scaled Jacobi coordinates defined as 




(7.1) 




(7.2) 
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and 



y=VJI^)[r. ^^j, (7.3) 

i.e., up to a scaling factor, the relative coordinate between the clusters 1 and 2 and the 
relative coordinate of their centre of mass with respect to the core. With Laplacians 
Ax and Ay with respect to x and y, the Hamiltonian ( 17.11 ) of this three-body projectile 
can be rewritten as 

Hq = -- — (A, + Ay) + y,i + y,2 + Vn- (7.4) 

To investigate the breakup cross sections for this system, we need wave functions at 
both positive and negative energies. 

In the notation of Refs. II103II105]| . the hyperradius p and hyperangle a are de- 
fined as 



p = y^x^+y^ (7.5) 

and 

a = arctan(y/;ii;). (7.6) 

The hyperangle a and the orientations £2^ and Qy of x and y provide a set of five 
angles collectively denoted as £2^. The volume element is dxdy = p^dpdQ^ with 
dQs = sin^ acos^ adadQ^dQy. 

The hyperspherical harmonics form an orthonormal basis which verifies a closure 
relation. The purely spatial hyperspherical harmonics read ||103N105ll 

<%S^5) = 4' ■ («) [YiA^.) ® Yi^Qy^]'''' . (7.7) 

where K is the hypermomentum quantum number, /_t and ly are the orbital quantum 
numbers associated with x and y, and L is the quantum number of total orbital mo- 
mentum. The functions 0^ ' depending on the hyperangle a are defined in Eqs. (9) 



and (10) of Ref. II105I . The hyperspherical harmonics involving spin are defined by 



JM 



(7.8) 



where Xs is a spinor corresponding to a total spin S of the three clusters. Intermediate 
couplings as, for example, the total spin si2 of the fragments are not displayed for 
simplicity. Index / stands for [IxlyLS). 

A partial wave function 0''*' is a solution of the Schrodinger equation ( 12.1b asso- 
ciated with the three-body Hamiltonian ( 17.41 ) at energy E. It can be expanded as 

</.^^(£,p,i25) = p-V2£;t^^(£,p) (7.9) 

yK 
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For bound states (E < 0), the hyperradial wave functions decrease asymptotically as 



x'yK iE,p)p-jM- yJ^'nN lEl/n'p). (7.10) 

Index T of Eq. ( 12. It is irrelevant for bound states within the present assumptions. 
The normalization of the scattering states (E > 0) is fixed by their asymptotic form. 
Several choices are possible. The asymptotic form of the hyperradial scattering wave 
function is for instance given by |fT9l 



where k ~ y Im^E /ti^ is the wave number and and are incoming and 
outgoing functions 01071 11081 11061 . In the neutral case, i.e. when clusters 1 and 
2 are neutrons, these functions read H^{x) ~ ±/(7rx/2)'''^ [7;^ (x) ± iT/f (jc)] where 
7/f and Yf^ are Bessel functions of first and second kind, respectively. In the charged 
case, expression ( 17.1 It is only an approximation because the asymptotic form of the 
Coulomb interaction is not diagonal in hyperspherical coordinates II1091 11061 . The 
indices Ja^a denote the partial entrance channels for this solution. For scattering 
states, index t of Eq. ( 12.1b is necessary and rather complicated: it represents the 
entrance channel ya^^a- The asymptotic behaviour of a given partial wave depends 
on the collision matrix. For real interactions, the collision matrix U'' of each partial 
wave J is unitary and symmetric. For three-body scattering, it differs from two-body 
collision matrices in an important aspect: its dimension is infinite since the particles 
can share the angular momentum in an infinity of ways. In practical calculations, 
its dimension depends on the number of hypermomenta included in the calculation, 
limited to a maximum K value, denoted as /iTmax- 

The three-body final scattering states are described asymptotically with two rel- 
ative wave vectors. Let k^, ki, k2 be the wave vectors of the core and fragments in 
the projectile rest frame. The asymptotic relative motions are defined by the relative 
wave vector of the neutrons 

k2i = VJiilk.v = (7.12) 

mi 

and the relative wave vector of the core with respect to the centre of mass of the 
fragments 

mi2kc. -»v(ki -f k2) 



kc(i2) = x/Mc(i2)k3, = — . (7.13) 



The total internal energy of the projectile with respect to the three-particle threshold 
is given by 
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E^-—k^^—{kl + kl). (7.14) 

The orientations of and 12^^, of k, and the ratio = arctan(yj^/ji[:j^) form the 
wave vector hyperangles 125^.. The hyperangle a<. controls the way the projectile 
energy E is shared among the fragments. For example, the energy of the relative 
motion between fragments 1 and 2 is £ cos^ . In the scattering states ( 12.2b . notation 
k^ thus represents Q^i^ and the final orientations Vc Vi, V2 of the three spins. It is 
convenient to replace these orientations by the total spin in of the fragments, the 
total spin S and its projection v. Relation ( I2.3l l is then given by 

X (^^^)^rwJE,p ,Qs). (7.15) 

where Ya = {lxa),ly(o,La),S). These functions are normalized with respect to 5(kx — 

kx)5(ky-ky)^SS'5vv'- 

The hyperradial wave functions XyK ^'^ determined from the Schrodinger 
equation ( 12. 11 1. The parity n = (—1)^ of the three-body relative motion restricts 
the sum over K to even or odd values. Rigorously, the summation over yK in ( I7.9l l 
should contain an infinite number of terms. In practice, this expansion is limited by 
the truncation value Km^x- The and values are limited hy lx + ly<K < Kmax- For 
weakly-bound and scattering states, it is well known that the convergence is rather 
slow and that large K^ax values must be used. 

The functions XyK derived from a set of coupled differential equations II105I 



rp- f cf- {K + 3/2){K + 5/2) 



2mN \dp^ p^ 



XyKiE,p) 



YK' 

where the potentials matrix elements are defined as 

yYK',rKiP)^i^yK'{^5)\ t Vijirj-r,)\?Vyf{Qs)). (7.17) 

i>j=l 

For bound states, approximate solutions can be obtained with an expansion on a fi- 
nite square-integrable basis. However, using such a basis for scattering states raises 
problems since they do not vanish at infinity. Their asymptotic form requires a 
proper treatment. This technical difficulty can be solved within the TJ-matrix theory 
111 101 11061 l55l which allows matching a variational function over a finite interval 
with the correct asymptotic solutions of the Schrodinger equation. 

In the TJ-matrix approach, both bound and scattering hyperradial wave func- 
tions are approximated over the internal region by an expansion on a set of square- 
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integrable variational functions defined over [Q,a]. Lagrange-mesh basis functions 
are quite efficient for describing two-body bound and scattering states |f56l II 1 II 
II 121 . The main advantage of this technique is to strongly simplify the calculation 
of matrix elements ( 17.17b without loss of accuracy if the Gauss approximation con- 
sistent with the mesh is used 11051 . This method was extended to three-body bound 
states in Ref. 01051 and to three-body scattering states in Ref. 01061 . We refer the 
reader to those references for details. 



7.2 Dipole strength distribution 

The El strength distribution for transitions from the ground state to the continuum 
is a property of the projectile that can be extracted from breakup experiments under 
some simplifying assumptions for cases where El is dominant. In the hyperspherical 
coordinate system, the multipole operators are given by Eq. (B2) of Ref. ill05l . For 
example, in two-neutron halo nuclei, the El strength is given by 

^-\p,n,)^eZ.'^P^Ym). (7.18) 

The El transition strength ( 13.161 ) from the ground state at negative energy £0 with 
total angular momentum Jq to the continuum is given by 



dB{El) _ 4 /mN 



5k 



dE 270 + 1 V fi^ 



(7.19) 



The El strength presents the advantage that it can also be calculated in various 
ways without constructing the complicated three-body scattering states 111 131 . Most 
model calculations of the El strength for ^He indicate a concentration of strength 
at low energies E lfT9ll2n[T07lll 14l[l 15111 161 . The origin of this low-energy bump 
remains unclear and can sometimes be attributed to a three-body resonance 01 071 
[191 . The existence of such a bump does not agree with the GSI data 111 171 . 

This puzzling problem deserves further studies. A first-order description of 
Coulomb breakup for ^He is probably not very accurate (see Sec. |7.3l l. even at 
the energies of the GSI experiment 111 17| . Extracting the El strength from breakup 
measurement is very difficult and not without ambiguities. This is exemplified by 
the variety of experimental results obtained for the breakup of the "Li two-neutron 
halo nucleus. As shown in Fig.[TT] most early experiments Illl8llll9lll20l did not 
display a significant strength at low energies in contradiction with data from the 



more recent RIKEN experiment 1121 
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Fig. 11 Experimental El strength for the breakup of the ' ' Li two-neutron halo nucleus: Ref. 11211 
(full circles), Ref. fTTSl (dash-dotted line), Ref fTT9l (histogram), Ref fT20l (dashed line). 
Reprinted figure with permission from Ref. (JJJJ. Copyright (2006) by the American Physical 
Society. 



7.3 The CCE approximation for three-body projectiles 

We consider a collision between a three-body projectile and a structureless target 
with mass mj and charge Zj-e 1 191 . The breakup reaction is described by the four- 
body Schrodinger equation ( 12.4b where /?o is given by Eq. (17.4b . The effective po- 
tential ( 12.5b between projectile and target is defined as 

y/.r(R,x,y) -y,r R+ — +Vxt R- — 



+V2T R 



nic y mi x \ 




In this expression, each interaction Vjp between a constituent of the projectile and 
the target is simulated by a complex optical potential (including a possible Coulomb 
interaction taking the cluster extension into account). 

In order to obtain breakup cross sections, one must calculate transition matrix 
elements for the breakup into three fragments. The transition matrix elements (14.11b 
read 



Tf, = {linlic(n))'^'^ 

A-) 



ralsv , P , ) I (R, X, y ) I f ) (R, p , i25 )) (7.21) 



for four-body breakup. The factor (Ati2Mc(i2)) appears when the integration is 
performed in coordinates p and Cl^ and the bound-state wave function ( 17.91 ) is 
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normed in this coordinate system rather than in Jacobi coordinates |fT9l . At the 
eikonal approximation, the exact scattering wave function f in Eq. (17. 21b is re- 
placed by its approximation given by Eqs. ( 14.11 ) and ( 14.171 ). The transition matrix 
element ( 17.211) is then obtained following Eq. ( 14.131 ) as 

Tf, = ihv j t/be-»>''4^''^(£,f25,,b). (7.22) 

with the eikonal breakup amplitude ( 14.141 ). that reads here |fT9l 

C(£,i25,,b) = (Mi2M.(i2))-'/' 

X {ip'^^l,^{E^p,£l,)\e'^^''^^-^^ (7.23) 

Following Eq. (17.20b . the eikonal phase shift % defined in Eq. ( 14.18b is obtained as 

X^XcT + X\T + X2T- (7.24) 

It depends on the transverse part b of R as well as on the transverse parts Sx and Sy 
of the scaled Jacobi coordinates x and y. 

From the transition matrix elements (17.21b . various cross sections can be derived. 
The differential cross section ( 14.15b with respect to the eight independent variables 
Q., k2i, kc(i2) reads in the cm. frame 

do ^ 1 ^_(}Ly^y\T\^ (125) 



The physical wave numbers k2\ and A:^(i2) are proportional to and ky and can thus 
be expressed from k and [ 19 |. Integrating Eq. ( 17.25b over all angles Q and Q.^^ 
leads to the energy distribution cross section da/dE. 

The CCE approximation has allowed calculating various elastic and breakup 
cross sections for ^He on ^''^Pb by treating *He as an a + n + n three-body sys- 
tem lfT9ll . In Fig. [12] the contribution from the different partial waves is displayed at 
240 MeV/nucleon. As expected for a transition from a 0+ ground state, the / = 1 ^ 
component is dominant. However the J — 0^ and 7 = 2+ components are not neg- 
ligible. The known 2+ resonance at 0.82 MeV is clearly visible in the total cross 
section. Extracting an El strength from such data is thus not easy, even at this high 
energy. 

A comparison of the CCE cross section (full line) with GSI data 11 171 is pre- 
sented in Fig. [13] The disagreement already discussed for the El strength in Sec. 17. 21 
is clearly visible. The data do not show as large a cross section at low energies as 
the theory. It is not even clear whether the 2+ resonance is visible in these data. 
Nevertheless the agreement is reasonably good above 2 MeV, given that no param- 
eter is fitted to this experiment in the model calculation. The 1^ contribution is 
calculated with two different ways of eliminating the unphysical bound states in the 
a -I- n potentials (dashed and dotted lines). The low-energy peak corresponds to a 
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Fig. 12 CCE calculation of the total and 0+, 1^, 2+ partial cross sections of ^He breakup on ^"^Pb 
at 240 MeV/nucleon |19| . Reprinted figure with permission from Ref. |19| . Copyright (2009) by 
the American Physical Society. 



broad resonance in the lowest 1 three-body phase shift. Further experimental and 
theoretical works are needed to explain this discrepancy. 
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Fig. 13 Comparison 1191 between the total CCE cross section (full line) of ^He breakup on ^"'^Pb at 
240 MeV/nucleon with the experimental data of Ref |117|. The 1^ partial cross sections calculated 
with two types of elimination of forbidden states (supersymmetry: dashed line, projection: dotted 
line) are also displayed. Reprinted figure with permission from Ref. (19]. Copyright (2009) by the 
American Physical Society. 

The advantage of the relative simplicity of the CCE is that various types of an- 
gular differential cross sections can be calculated. Examples of double differential 
cross sections showing various energy repartitions between the fragments are pre- 
sented in Fig. 7 of Ref. IH. 
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7.4 The CDCC approximation for three-hody projectiles 

The CDCC method has also been extended to three-body projectiles. In the first 
works 1521 l60l l53l . the pseudostate discretization was adopted. Indeed, it avoids 
the difficult construction of scattering states and allows an accurate treatment using 
expansions involving Gaussians with various widths. Only recently was the con- 
struction of bins attempted ||2T1 . The difficulty of the calculation restricted the first 
applications to elastic scattering. 

The differential cross section for elastic scattering of ^He on '"C at 229.8 MeV 
is displayed in Fig. [141 A single-channel calculation neglecting breakup channels 
(dotted line) overestimates the experimental data of Ref. 11221 . The shape of the 
data is very well reproduced by introducing 0+ and 2+ pseudochannels and taking 
account of all couplings (full line). In Fig. [151 is displayed a comparison between 
calculations of ^He elastic scattering on ^"^Bi at 22.5 MeV involving two-cluster 
(three-body, dashed line) and three-cluster (four-body, full line) descriptions of ^He. 
A significant difference appears between calculations neglecting breakup channels 
("no coupling") and those including it ("full coupling"). The agreement with exper- 
imental data 11231 11241 seems better within the four-body treatment including the 
breakup channels. 




Fig. 14 Ratios of differential cross sections obtained with CDCC to Rutherford cross section for 
the elastic scattering of 'He on '^C at 229.8 MeV without and with coupling to breakup chan- 
nels 1521 . Experimental data from Ref 11221 . Reprinted figure with permission from Ref. (52). 
Copyright (2004) by the American Physical Society. 

Another type of basis functions, based on deformed oscillators, has been used to 
construct ^He pseudostates in Ref. Il53l . This technique also allowed a description 
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Fig. 15 Ratios of differential cross sections to Rutherford cross section for the elastic scattering of 
^He on2*''^Biat22.5 MeV: comparison of three- and four-body CDCC without and with coupling to 
breakup channels 1601 . Experimental data from Refs. I123II1241 . Reprinted figure with permission 
from Ref. 1601 . Copyright (2006) by the American Physical Society. 



of elastic scattering explicitly including breakup channels. In Fig. [16] the elastic 
scattering of ^He on ^^Zn at 13.6 MeV is compared with experimental data from 
Ref. II125I . These results show that including partial waves up to 7 = 2 and taking 
coupling into account (full line) allow a good agreement with data (dots). Here also, 
the calculation omitting the coupling to the continuum (dashed line) disagrees with 
the experimental data. The same basis has recently been extended to the construction 
of bins ED. 

While, for three-body projectiles, the effect of breakup channels has been in- 
cluded for some time in studies of elastic scattering, the determination of breakup 
cross sections is just starting. Some preliminary calculations have been published 
recently. Some of them are not fully converged |!2TI| or involve simplifying assump- 
tions II20I . A recent CDCC calculation 1611 provides a good agreement with experi- 
ment HI 171 for ^He breakup on '"C. For ^He breakup on ^''^Pb, it does not agree well 
with experiment and is about a factor of two lower than the CCE results of Ref. fTg) 
displayed in Fig.[T3] The reasons of these discrepancies are not yet understood. Nev- 
ertheless, the CDCC method should allow a precise treatment of three-body breakup 
in a near future. 
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Fig. 16 Ratios of differential cross sections to Rutherford cross section for the elastic scattering 
of '■He on '''^Zn at 13.6 MeV: comparison of CDCC calculations with various basis sizes and max- 
imum energies £|„ax with a single-channel calculation (55] ■ Experimental data from Ref. [125]. 
Reprinted figure with permission from Ref. 1531 . Copyright (2008) by the American Physical So- 
ciety. 



8 Perspectives 



The theory of breakup reactions offers several accurate non-relativistic approxima- 
tions covering a broad energy range, that allow an interpretation of various exper- 
iments. A good accuracy is reached for some time for the breakup of two-body 
projectiles and is in view for the breakup of three-body projectiles. Good results 
can already be obtained with the simplest models of projectile structure, provided 
that the value of the projectile binding energy is correct. This suggests that only 
limited spectroscopic information can be extracted from the comparison of theory 
and experiment. This is partly due to the fact that a comparison of experimental 
data with results of calculations usually requires complicated convolutions. Nev- 
ertheless, breakup has proved to be an efficient alternative probe to measure the 
separation energies of bound states of exotic nuclei. When performed on light tar- 
gets, it also provides information about the location and width of resonances of 
such nuclei. Moreover, some information about the quantum numbers of the ground 
state of exotic nuclei can be assessed from breakup measurements. The extraction 
of spectroscopic factors, however, is very sensitive to the accuracy of the absolute 
normalization of experiments. Moreover, the sensitivity of breakup calculations to 
the description of the continuum of the projectile indicates that these extractions 
should be performed with caution. In addition. Coulomb breakup on heavy targets 
is also used to measure astrophysical S factors. However the accuracy of this indirect 
technique is uncertain. 
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Several methods can now be applied to the breakup of three-cluster projectiles 
(CDCC method, eikonal approximation, . . . ). They will allow studying coincidence 
observables that are more difficult to measure but less sensitive to the absolute nor- 
malization of cross sections. They should also allow the study of correlations be- 
tween the emitted fragments. In this respect, efforts should be made at the inter- 
face between theory and experiment to facilitate the transformation of the results 
of model calculations into quantities comparable with the data, taking account of 
the resolution and acceptances of the detection setup. On the theoretical side, three- 
cluster bound states can be obtained with good accuracy but the difficult treatment 
of the three-body continuum still requires progress. 

Attempts to improve the model description of the projectile by including ex- 
cited states of the clusters composing the projectile have started with the extended 
CDCC. In the future, one can expect a further improvement by using a microscopic 
description of the projectile within the microscopic cluster model II1261 11271 1128L 
involving effective nucleon-nucleon forces and full antisymmetrization. Improve- 
ments in the projectile description should first concern bound states. This should 
reduce the uncertainties appearing in non-microscopic cluster models because of 
the effective forces between the clusters in the projectile and between the clusters 
and the target. Using fully antisymmetrized wave functions in breakup calculations 
seems to be within reach for two-cluster projectiles. This approach should open the 
way towards ab initio descriptions of the projectile based on fully realistic nucleon- 
nucleon forces. 

All the reaction descriptions presented in this review have been developed within 
non-relativistic quantum mechanics. However, relativistic effects may be significant 
and affect the analysis of breakup data, even at intermediate energies of a few tens 
of MeV/nucleon. Several authors have started analyzing these effects and have pro- 
posed ways to take them into account in time-dependent 11291 or CDCC II130II13I1 
frameworks. Since some of the new facilities of radioactive-ion beams will operate 
at high energies (a few hundreds of MeV/nucleon), these effects will have to be 
better understood and incorporated in state-of-the-art reaction models. 
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